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ABSTRACT 

The mysterious “21 /xrri” emission feature seen almost exclusively in the short-lived 
protoplanetary nebula (PPN) phase of stellar evolution remains unidentified since its 
discovery two decades ago. This feature is always accompanied by the equally mys¬ 
terious, unidentified “30 /xm” feature and the so-called “unidentified infrared” (UIR) 
features at 3.3, 6.2, 7.7, 8.6, and 11.3 /xm which are generally attributed to polycyclic 
aromatic hydrocarbon (PAH) molecules. The 30 /xrn feature is commonly observed in 
all stages of stellar evolution from the asymptotic giant branch (AGB) through PPN to 
the planetary nebula phase. We explore the interrelations among the mysterious 21 /xm, 

30 /xm, and UIR features of the 21 /xrn sources. We derive the fluxes emitted in the ob¬ 
served UIR, 21 /xm, and 30 /xm features from published ISO or Spitzer/ IRS spectra. We 
find that none of these spectral features correlate with each other. This argues against a 
common carrier (e.g., thiourea) for both the 21 /xm feature and the 30/xm feature. This 
also does not support large PAH clusters as a possible carrier for the 21 /xm feature. 

Subject headings: circumstellar matter — dust, extinction — infrared: stars — stars: 
AGB and Post-AGB — stars: evolution 


1. Introduction 

In carbon-rich evolved objects, there are two prominent, mysterious emission bands known as 
the “21/xm” and “30/xm” spectral features (see Jiang et al. 2010). The enigmatic 21 /xm emission 
feature was first detected in four protoplanetary nebulae (PPNe) through the 7.7-22.6 /xm spectra 
obtained by the Low Resolution Spectrometer (LRS) on board the Infrared Astronomical Satellite 
(IRAS; Kwok et al. 1989). To date, this feature has been seen unambiguously in ~ 18 objects in 
the Milky Way (Cerrigone et al. 2011) and ~9 objects in the Large and Small Magellanic Clouds 


department of Physics and Astronomy, University of Missouri, Columbia, MO 65211, USA; 

amishra@mail.missouri.edu, lia@missouri.edu 

department of Astronomy, Beijing Normal University, Beijing 100875, China; bjiang@bnu.edu.cn 



2 


(Volk et al. 2011)0 The 21 /im sources are all PPNe, objects in a rapid transition period (of several 
thousand years) between the asymptotic giant branch (AGB) and planetary nebula (PN) phases0 
With a peak wavelength at ~20.1/un and a FWHM of ~ 2.2-2.3/mi, the 21/mi feature displays 
little shape variation among different sources. This feature emits up to ~8% of the total infrared 
(IR) power of a 21 /mi source. The carrier of this feature remains unidentified, despite the fact that 
over a dozen carrier candidates have been proposed, including inorganic materials: TiC (von Helden 
et al. 2000, but see Li 2003), SiS 2 (Goebel 1993), SiC (Speck & Hofmeister 2004, but see Jiang 
et al. 2005), Si02-coated SiC (Posch et al. 2004), carbon-silicon mixtures (Kimura et al. 2005), 
FeO (Posch et al. 2004, but see Li et al. 2013), Fe 203 and Fe 304 (Cox 1990, but see Zhang et al. 
2009a), and organic materials: Ti-coordinated caged carbon (Kimura et al. 2005), urea or thiourea 
(Sourisseau et al. 1992), polycyclic aromatic hydrocarbon (PAH), and hydrogenated amorphous 
carbon (HAC) (Buss et al. 1993, Justtanont et al. 1996). It is considered as “one of the most 
interesting unresolved mysteries in astrochemistry ” (Kwok et al. 2002). 


The 30/mi feature, first discovered by Forrest et al. (1981) in several carbon stars and in two 
PNe through the Kuiper Airborne Observatory spectrometry, is very broad and strong, extending 
from ~ 24 /mi to ~ 45 /mi and accounting for up to ~ 30% of the total IR luminosity of a 30 /mi 
source (Volk et al. 2002). Unlike the 21 /mi feature which displays little shape variation, the 30/mi 
feature varies in its peak wavelength and width among different sources (e.g., see Hrivnak et al. 
2000, Hony et al. 2002). 


The carrier of the 30 /mi feature also remains unidentified. Magnesium sulfide (MgS) solids 
have long been proposed as a candidate carrier (Goebel & Moseley 1985; Nuth et al. 1985; Jiang 
et al. 1999; Szczerba et al. 1999; Hony et al. 2002, 2003; Lombaert et al. 2012). However, Zhang et 
al. (2009b) recently ruled out MgS as a valid carrier since it would require too much S to account 
for the observed large amount of power emitted from the 30 /mi feature (also see Messenger et al. 
2013, Otsuka et al. 2014)0 


While so far the 21 /mi feature appears to be restricted to PPNe, the 30 /mi feature is more 
ubiquitously seen in a large number of carbon-rich objects at various evolutionary stages, including 
AGB stars, post-AGB stars and PNe (see Jiang et al. 2010, Zhang &; Jiang 2008). We note that the 
30 /mi feature is seen in all 21 /mi sources, although not all 30 /mi sources emit at the 21 /im feature. 
Furthermore, all 21 /im sources display the so-called “unidentified IR” (UIR) features at 3.3, 6.2, 


1 In the Magellanic Clouds the total number of 21 /an sources would be 13 if the marginal detection in the other 
four sources is included (Volk et al. 2011). 

2 A few weak detections of the 21 /im feature have also been reported for their precursors - highly evolved carbon 
stars (Volk et al. 2000) and their successors - PNe (Hony et al. 2001, Volk 2003). These detections are yet not very 
conclusive as the 21 /im feature claimed for these sources is considerably weak. 

3 Lombaert et al. (2012) modeled the 30 fim feature of the extreme carbon star LL Peg with MgS dust. They 
argued that if MgS is in thermal contact with amorphous carbon and SiC, the amount of MgS required to reproduce 
the strength of the 30 /im feature agrees with the solar abundance of sulfur. 
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7.7, 8.6 and 11.3 /im (Hrivnak et al. 2008) which are commonly attributed to PAH molecules (Leger 
& Puget 1984, Allamandola et al. 1985). 

To gain insight into the nature of the carriers of the 21 fim and 30 /im features, we select ten 
well-studied Galactic 21 /im sources. The 21 /im, 30 /im and UIR features of these sources have 
been observationally obtained. This allows us to derive the fractional fluxes emitted respectively 
in the 21 /im, 30 /im and UIR features (see f|2]). We examine in the interrelations among the 
21 /im, 30 /mi and UIR features in the selected 21 /im sources. We discuss in )J5]the implications of 
these interrelations. Ij5] discusses the methods of deriving the fluxes of the 21 /im, 30/im and UIR 
features. We summarize the major results in Sj6l 


2. Deriving the Fluxes of the 21 /im, 30 /im and UIR Features 

We have selected ten well-studied Galactic 21 /im sources. For these sources, the high-quality 
mid-IR spectra obtained by the Spitzer Space Telescope and the Infrared Space Observatories (ISO) 
allow one to measure relatively accurately the (integrated) fluxes emitted in the 21 /im feature (F 21 ), 
in the 30 /mi feature (F 30 ), and in the UIR features (Fuir)EI In Table CD we tabulate the stellar 
and circumstellar parameters of the ten Galactic 21 /im sources (compiled from the literature): the 
effective temperature T e g , luminosity L*, core mass M*, stellar radius r*, and distance d from Earth 
of the central star. 

Below we briefly comment on the individual sources, focusing on the nebular morphology of 
each source. IRAS 02229+6208 is a cool, highly reddened post-AGB star. It has an elliptically 
extended nebula as revealed by the polarization map of Ueta et al. (2005). IRAS 04296+3429 is in 
the advanced post-AGB evolution stage with a bipolar lobe structure as revealed by the scattered 
light images at 0.56 and 0.81 /im (Sahai 1999). IRAS 05341+0852 has a very extended, optically 
thin atmosphere. Its visible image shows an elongated elliptical nebula around the central star 
(Ueta et al. 2000). IRAS 07134+1005 (HD 56126) is one of the best studied post-AGB star. It 
has an axial symmetric structure (Meixner et al. 1997; Ueta et al. 2000). It is one of the sources in 
which the 21 /im and 30/trn features were first discovered (Forrest et al. 1981, Kwok et al. 1989). 
IRAS 16594-4656 is a bipolar post-AGB star as indicated by its dust emission spectral energy 
distribution (Meixner et al. 1999; Hrivnak et al. 2008). It has an optically thick circumstellar 
envelope. IRAS 20000+3239 has an extended and axisymmetric bipolar structure as revealed by 
the near-IR imaging polarimetry of this object (Gledhill et al. 2001). For IRAS 22223+4327, the 
polarization map reveals an extended, optically thick circumstellar envelope (Gledhill et al. 2001). 
IRAS 22272+5435 is extremely carbon-rich and fairly bright both in the IR and optical (Ueta et 
al. 2001). IRAS 22574+6609 is very faint in the visible. It shows a bipolar morphology with a 


4 To have a full count of the 21 /im, 30 /tm, and UIR fluxes, we require that the IR spectra of the selected sources 
at least cover the wavelength range of ~ 5-35 /tm. As a result of this criterion, the 21 /tm sources which lack the IR 
spectra at A < 10 /tm are not included in this study. 
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dark lane dividing the nebula into two lobes (Ueta et al. 2000; Su et al. 2001). IRAS 23304+6147 
has a quadrupolar shape, although its two pairs of lobes are not well separated. This object is 
remarkably similar to IRAS 20000+3239 in its post-AGB properties and AGB ejecta (Sahai et al. 
2007). 

We retrieve the IR spectra of the selected sources from the Spitzer and ISO archives. All 
sources exhibit the 21 pm, 30 pm and UIR features. In order to quantify the relative contributions 
from different components, we perform spectral decomposition of the IR spectra of these sources. 
Following Zhang et al. (2010), we use the IDL package PAHFIT of Smith et al. (2007) to decompose 
the observed IR spectra into individual components, including contributions from stellar continuum, 
H 2 lines, thermal dust continuum, and the 21 pm, 30 pm and UIE features. This decomposition 
technique will be discussed in (jH We note that the decomposition is not really a physical dust 
model. It only serves to measure the fluxes in the features. 

The stellar continuum is assumed to be a blackbody of T e g, the stellar effective temperature. 
In all cases, the stellar contribution to the IR spectra is negligible. Following Zhang et al. (2010), we 
include eight H 2 S(0)-S(7) rotational lines because they are seen in the observed spectra and need to 
be subtracted in order to accurately determine the fluxes of the dust spectral features. The profiles 
of the H 2 lines are assumed to be Gaussian. A modified blackbody I\ A 2 B\{T) is adopted 
to simulate the dust thermal continuum, where B\(T) is the Planck function of temperature T at 
wavelength A. For all the sources the best-fit requires two components with different temperatures 
to describe the dust thermal continuum: a warm dust component with T\y ~ 130-190 K, and a 
cold component with Tq ^ 60-90 K. The 21 pm, 30 pm and UIR features are fitted with a number 
of Drude profiles I\ = Io7 2 / {(A/Ao — Ao/A) 2 + 7 2 }, where Ao is the central wavelength, Iq is the 
central intensity, and 7 is the fractional FWHM of each feature. The broad 30 pm feature is fitted 
by two or three subfeatures at 26, 30 and 33 pm§ For the UIR features, we include all the PAH 
bands listed in Table 1 of Draine &; Li (2007). The Drude profile is expected for classical damped 
harmonic oscillators (see Li 2009). It closely resembles a Lorentzian profile. Compared with a 
Gaussian profile, both the Drude profile and the Lorentzian profile have more extended wings. The 
21 pm, 30 pm and UIR features are better fitted with Drude profiles than Gaussian profiles, while 
the H 2 lines are better fitted with Gaussian profiles. 

The optimal fitting to the observed spectra is achieved through the Levenberg-Marquardt least- 
square algorithm implemented in the IDL package. We are able to obtain reasonably good fits to the 
observed spectra using this technique. In Figures [I] [3] we show the spectral decomposition fitting 
to the ten sources. The fluxes emitted by the 21 pm, 30 pm and UIR features are listed in Table [2j 


5 The ISO/SWS spectroscopy suggested that the 30 pm feature consists of two subfeatures: a narrow feature at 
26 pm and a broad one at 33 pm (Hrivnak et al. 2000, Volk et al. 2002). However, the two distinct components in the 
30 pm feature are less prominent in the Spitzer spectra (see Hrivnak et al. 2009). Like Zhang et al. (2010), we find 
it practical to use two or three Drude profiles peaking at 26, 30 and 30 pm to represent this feature as they provide 
better fits to the data. 
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Also listed in Table [2] is -Fir, the dust continuum IR flux integrated over the 5-45 /im wavelength 
ranged The uncertainties of the feature strengths are estimated using the full covariance matrix of 
the least square parameters as described in Smith et al. (2007). 


3. Interrelations among the 21 /xm, 30 /mi and UIR Features 

Let M 2 i urn , A/ao^rn, and Muir respectively be the mass of the carrier of the 21 /im, 30/xm, and 
UIR features. Let and respectively be the ultraviolet/visible mass absorption 

coefficient of the carrier of the 21 /im, 30 /im, and UIR features. It is reasonable to assume that 
the total masses of the carriers of these unidentified features (M2i /im , . and Mrir) are 

proportional to the total mass of the bulk dust (M dust ) in the shell, although the proportional factor 
may differ from one to another. Let / 21 , /30, and /uir respectively be the ratio of the mass of the 
carrier of the 21 /tm, 30 /xm, and UIR features to the mass of the bulk dust in a 21 gm source. With 
Fir oc (Mdust x F*) we have F 2 \ oc ^M 21/im x K^ m x L *) (/21 x Af dust x F*) oc (/ 21 x Fir). 

Similarly, we have F30 oc (/30 x Fir) and Frir oc (/uir x Fir). Therefore, even if /30 does not 
correlate with /21 or /uir, F30 could still correlate with F 2 \ and Fjir through Fir. Therefore, 
for a physically reasonable correlation study of the 21 /an, 30 /im and UIR features, one need to 
normalize F 2 1, F30 and Frir by Fir, in order to cancel out their common proportionality to Fir. 

We explore the interrelations among the 21 /xm, 30/xm and UIR features. We normalize the 
fluxes emitted from these features (F 2 i , F30, Frir) by Fir, the total near- to mid-IR emission in the 
wavelength range of 5-45 /xm. This is to see what fractions of the total near- to mid-IR luminosity 
are being carried by these features. 

We show in Figure |4] the interrelations of the fractional fluxes of the 21/xm, 30 /xm and UIR 
features. It is clear from Figure[J]that these features are not correlated. This differs from Cerrigone 
et al. (2011) who reported a positive correlation between Frir and F 2 i and suggested that the 
abundances of the carrier of the UIR bands and of the carrier of the 21 /xm feature may be linked 
together, possibly because of a common formation process. We attribute the difference to the fact 
that Cerrigone et al. (2011) did not normalize Frir and F 2 1 by Fir. We argue that the correlation 
between F 2 i and Frir found by Cerrigone et al. (2011) may merely reflect the fact that both F 2 i 
and Frir correlate with F m (i.e., with M 21fim oc M dust , M 30fnri oc M dust , and M mR oc M dust , even 
if /30 does not correlate with f 2 \ or /uir, could still correlate with M 2 i^ m and Mrir). 


®For IRAS 04296+3429, IRAS 05341+0852, and IRAS 22223+4327, the IR spectra are not available at A > 35 /xm. 
We derive Fir for these sources by extrapolating the observed spectra to A = 45 /xm. 
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Table 1: Stellar Parameters for the Ten Selected 21 (tm Sources. 


IRAS 

Sources 

F e ff 

(K) 

F* 

(F©) 

A/* 

(AF 0 ) 

T* 

(Re) 

d 

(kpc) 

02229+6208 

5500 

8333 

0.558 

98.4 

2.1 

04296+3429 

6500 

8333 

0.554 

69.1 

5.0 

05341+0852 

6500 

8430 

0.551 

69.5 

7.8 

07134+1005 

7250 

6554 

0.841 

49.0 

2.2 

16594-4656 

10000 

10279 

0.60 

32.0 

2.55 

20000+3239 

5500 

5186 


77.0 

2.24 

22223+4327 

6500 

6075 

0.551 

59.0 

3.2 

22272+5435 

5650 

10990 

0.574 

106.0 

1.55 

22574+6609 

5500 

8061 

0.60 

96.0 

7.6 

23304+6147 

6750 

8347 

0.66 

64.0 

3.25 


Table 2: Integrated Fluxes Emitted in the UIR (Fuir), 21 (tm (F 21 ), and 30/im (F 30 ) Features as well as the Total 
Near- to Mid-IR Emission (Fir) in the ~ 5-45 (tm Wavelength Range, as Derived from the PAHFIT Decomposition 
Method._ 


IRAS 

Sources 

Warm Dust 
T w (K) 

Cold Dust 
Tc (K) 

UIR 

(10 ~ 12 Wm" 2 ) 

F21 

(10 ~ 12 Wm“ 2 ) 

F30 

(10 ~ 12 Wm~ 2 ) 

Fir 

(10“ 12 Wm- 2 ) 

02229+6208 

160 

90 

5.78+0.016 

0.322+0.002 

9.10+0.138 

25.6+1.95 

04296+3429 

185 

61 

0.870+0.006 

0.302+0.002 

1.62+0.012 

7.15+0.23 

05341+0852 

190 

70 

0.283+0.005 

0.023+0.008 

0.35+0.002 

1.66+0.07 

07134+1005 

160 

90 

1.84+0.017 

1.40+0.002 

2.88+0.015 

17.3+1.60 

16594-4656 

135 

78 

3.25+0.014 

2.57+0.002 

9.97+0.080 

43.6+2.90 

20000+3239 

158 

80 

0.522+0.009 

0.185+0.001 

2.63+0.019 

7.29+1.50 

22223+4327 

150 

70 

0.158+0.005 

0.114+0.001 

0.986+0.014 

5.42+0.30 

22272+5435 

165 

90 

4.22+0.036 

0.840+0.010 

11.4+0.135 

35.8+2.20 

22574+6609 

170 

60 

0.734+0.008 

0.144+0.002 

1.20+0.006 

4.61+0.69 

23304+6147 

165 

89 

0.767+0.011 

0.416+0.001 

2 . 88 + 0.010 

8.90+0.83 


Note. — For IRAS 04296+3429, IRAS 05341+0852, and IRAS 22223+4327, the IR spectra are not available at 
A > 35 (tm. We derive Fir for these sources by extrapolating the observed spectra to A = 45 (tm. 



7 





i 111 1 1 111 1 1 111 1 1 111 1 1 111 1 1 

15 20 25 30 35 40 

A +m) 



Fig. 1. — Decomposition of the observed IR spectra (black points) of the 21 /im sources (IRAS 04296+3429, 
IRAS07134+1005, IRAS20000+3239, and IRAS22272+5435) into a stellar continuum (pink line), two dust ther¬ 
mal continuum emission components (red lines), the 21 /im, 30/im and UIR features (blue lines), and the Ho bands 
(violet lines). The decomposition is done with the PAHFIT technique. Gray lines are the sum of the stellar and dust 
thermal continuum. Green lines show the resulting model spectra. 
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Same as Figure[l]but for IRAS 02229+6208, IRAS 05341+0852, IRAS 16594-4656, and IRAS 22223+4327. 



Fig. 3.— Same as Figure[l]but for IRAS22574+6609 and IRAS23304+6147. 
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Fig. 4.— Interrelations among the 21 pm, 30 /im, and UIR features of ten Galactic 21 pm sources: (a) F 30 vs. F 21 , 
(b) -F 30 vs. Fuir, and (c) F 21 vs. Fuir. All quantities are normalized by Fir, the total near- to mid-IR emission 
obtained by ISO or Spitzer in the 5-45 pm wavelength range, to cancel out their common proportionality to Fir. The 
fluxes emitted in the 21 pm, 30 pm, and UIR features are derived from the PAHFIT decomposition method. 
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4. Implications 

The UIR bands are commonly attributed to PAH molecules. If the 21 gn i feature is due to 
large PAH clusters (Buss et al. 1990), F 21 and T\jir should be somewhat correlated: the fast stellar 
wind from the central star during the PPN phase might shatter PAH clusters into PAH molecules, 
possibly obliterating the carrier of the 21 gm feature. The lack of correlation between Turn and 
T 21 indicates that the 21 (im feature is probably not due to large PAH clusters. 

Papoular (2000, 2011) hypothesized that macromolecules with thiourea functional groups and 
aliphatic chains (made of CH 2 groups, oxygen bridges and OH groups) attached to various car¬ 
bonaceous structures (mainly compact and linear aromatic clusters) are responsible for the 21 grn, 
30 gm, and UIR features, with thiourea emitting the 21 gm feature, the aliphatic chains producing 
the 30 gm feature, and the aromatic clusters giving rise to the UIR bands. If this is the case, 
one would expect the 21 /im, 30 /im and UIR features to be positively correlated. The lack of any 
positive correlation among these features as demonstrated in §2 does not support the hypothesis of 
Papoular (2000, 2011). 

Volk et al. (2011) reported an anti-correlation between the 30/im and UIR features for the 
Magellanic Cloud 21 gm sources. They postulated there may be some competition in forming the 
UIR and 30 gm feature carriers. Indeed, if the 30 gm feature is due to hydrogenated amorphous 
carbon (HAC), such an anti-correlation could result from radiation-induced decomposition of HAC 
grains into PAHs (Scott et al. 1997) and/or stellar wind-induced shattering of HAC grains into 
PAHs. However, as shown in Figure |4] for the Galactic 21 /im sources the 30 gm and UIR features 
do not seem to be anti-correlated. 

Spectroscopically, the Magellanic Cloud 21 gm sources differ appreciably from the Galactic 
21 gm sources (see Volk et al. 2011): (1) the average 21 gm feature strength for the Magellanic 
Cloud objects is smaller than that observed in the Galaxy; (2) the 30 gm feature of the Magellanic 
Cloud sources carries less of the total rnid-IR emission than is normally the case for the Galactic 
21 (im sources; (3) the UIR features of the Magellanic Cloud 21 gm sources are typically twice 
as strong as in the Galactic sources, with the peak of the UIR emission often being higher than 
that of the 30 /im feature; and (4) the Magellanic Cloud objects have more structure in their 
UIR spectra and show more variation in the shapes and strengths of the 7.7 and 8.6 gm features 
than is observed in the Galactic 21 /im sources; more particularly, most of the Galactic 21 gn 1 
sources have “unusual” UIR properties compared to the typical interstellar UIR bands (e.g., the 
Galactic 21 gm sources have a broad 8 (im feature rather than separately resolved 7.7/8.6 (im UIR 
features) while the Magellanic Cloud 21 (in 1 sources show more “normal” looking UIR features (see 
Peeters et al. 2002). The remarkable difference in the spectral appearance of the UIR features 
between the Galactic and Magellanic Cloud 21 gm sources would potentially hold important clues 
for understanding how the physical and chemical conditions affect the chemical structures of the 
carriers of these unidentified features. 

An unique characteristic of the mysterious 21 gm feature is that it is almost exclusively seen in 
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the short-lived PPN phase of stellar evolution and always accompanied by the equally mysterious 
30 /xm feature and UIR bands, while the 30 /xm feature is commonly observed in all late stages of 
stellar evolution from the AGB through the PPNe to the PNe phase (the UIR features are seen 
in PPNe and PNe, but not in AGB stars). This implies that the carriers of the 21 /xm feature are 
likely formed during the so-called superwind phase (e.g., see von Helden et al. 2000), which is a 
phase of high mass-loss where AGB stars lose the remaining stellar envelope, terminating their life 
on the AGB (see Renzini 1981). They may also be easily destroyed by the highly energetic photons 
available in PNe. We also note that while the UIR features are commonly seen in PPNe, PNe, and 
the ISM, they are usually not seen in the mid-IR spectra of carbon-rich AGB stars!! 

Finally, we also note that, in addition to the variable 30 /im feature (in contrast to the invariable 
21 /im feature)]! the overall shape of the Spitzer spectra of the 21 /xrn sources show also a lot of 
diversity, from having almost no flux at ~ 10 /xm (like IRAS 22223-4327) to having a very flat 
spectrum (like IRAS 05341+0852, IRAS 22574+6609). We do not know what causes this diversity. 
We have examined the stellar luminosity, and the optical depth, the expansion history, and the 
inner and outer edges of the dust shell. They do not seem to be responsible for this overall spectral 
shape diversity. 


5. PAHFIT Decomposition Versus Spline Fitting 

When we use the PAHFIT decomposition method to derive the fluxes emitted from the 21 /xm, 
30 /xm, and UIR features (see )J2]), we approximate the dust continuum with two modified blackbod- 
ies (see Figures 11121) . We admit that in dust spectral analysis the dust continuum determination is 
often notoriously difficult. This is true for studying PAHs (e.g., see Uchida et al. 2000), silicates 
(e.g., see Xie et al. 2014, Shao et al. 2014), and many other dust species. For the 21/xm sources 
discussed in this work, what is the true dust continuum? Apparently, the fluxes of the dust spectral 
features are sensitive to the continuum determination. The key question here is to what extent the 
assumed continuum will affect the interrelations among these features examined in )}3]? To examine 
this, instead of approximating the dust continuum as two modified blackbodies, we take several 
points at ~ 6 /xm, ~ 10 /xm, ~ 18 /xm, ~ 24 /xm, and ~ 44 /xm to define an underlying continuum and 
fit the continuum with a cubic spline function (see Figures I5I6D . The 21 /xm, 30/xm, and UIR fea- 


7 The few carbon stars that display the UIR features all have a hot companion that emits UV photons (Speck & 
Barlow 1997, Boersma et al. 2006). It has been argued that the carriers of the UIR features (e.g., PAHs) are present 
in carbon stars but are simply not excited sufficiently to emit at mid-IR due to lack of UV photons. However, Li & 
Draine (2002) have shown that the excitation of PAHs does not require UV photons, visible and near-IR photons are 
capable of exciting PAHs to emit at the UIR bands. The visible/near-IR cross sections measured by Mattioda et al. 
(2005) for PAH ions imply that PAHs can be excited by the soft stellar photons from C stars. Therefore, the absence 
of the UIR bands in C stars places a limit on the abundance of small PAHs in these outflows. 

8 Hony et al. (2002) found that the peak wavelength of the 30 /xm feature shifts from ~ 26 /xm in some AGB stars 
to ~ 38 /xm in PNs. 



^-F (10 12 W/m 2 ) AF (10 -'! w/m 2 ) 
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Table 3: Same as Table [2] but with the Dust Continuum Fitted with a Spline Function (See (J5l and Figures [51611 . 


IRAS 

Sources 

UIR 

(10 -12 W m -2 ) 

+21 

(10 12 Wm“ 2 ) 

+30 

(10~ 12 Wm" 2 ) 

+IR 

(10" 12 Wm“ 2 ) 

02229+6208 

5.05+0.57 

0.34+0.02 

7.00+0.58 

25.6+1.95 

04296+3429 

0.87+0.06 

0.35+0.02 

0.63+0.06 

7.15+0.23 

05341+0852 

0.28+0.04 

0.029+0.002 

0.15+0.01 

1.66+0.07 

07134+1005 

1.90+0.30 

1.62+0.06 

2.07+0.25 

17.3+1.60 

16594-4656 

3.46+0.41 

2.02+0.041 

4.10+0.07 

43.6+2.90 
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0.54+0.045 
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Fig. 5. — Same as Figure |T] but with the dust continuum fitted with a spline function instead of two modified 
blackbodies. 
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Fig. 6.— Same as Figure [2] but with the dust continuum fitted with a spline function instead of two modified 
blackbodies. 




Fig. 7. — Same as Figure [3] but with the dust continuum fitted with a spline function instead of two modified 
blackbodies. 
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tures are again fitted with Drude profiles. In Table [3] we list the fluxes emitted in the 21 /mi, 30 /im 
and UIR features derived from the spline continuum method, as well as the dust IR flux integrated 
over the 5-45 /mi wavelength range. They differ from that of the PAHFIT decomposition method by 
~ 10% for the 21 /mi and UIR features, and by ~ 30-60% for the 30 /mi feature. Figure [8] shows the 
interrelations of the fluxes of the 21 /mi, 30 /mi and UIR features (normalized by the 5-45 /mi dust 
emission) derived from the spline continuum method. Similar to that obtained from the PAHFIT 
decomposition method, we see no correlation among these features. Therefore, we conclude that 
although the true dust continuum is not precisely known and the true fluxes emitted in the 21 /im, 
30 /mi and UIR features are sensitive to the assumed continuum underneath the features, these fea¬ 
tures are unlikely correlated as demonstrated by the two very different continuum-determination 
methods. 


6. Summary 

We explore the interrelations among the mysterious 21 /im, 30 /mi and UIR features of the 
21 /mi sources. The principal results are the following: 

1. The 21/mi feature does not correlate with the 30 /mi feature. This argues against the hy¬ 
pothesis of thiourea and aliphatic chains attached to various carbonaceous structures as the 
carriers for both the 21 /mi feature and the 30 /mi feature. 

2. The 21 n m feature does not correlate with the UIR features. This argues against large PAH 
clusters as a possible carrier for the 21 /mi feature. 

3. The 30 /im feature and the UIR features are not correlated. This does not support the spec¬ 
ulation that the UIR carriers (e.g., PAHs) may result from the decomposition or shattering 
of the 30/mi-feature carrier (e.g., HAC). 


We thank Z. Ivezic, T. Ueta, K. Volk, Y.X. Xie, and K. Zhang for helpful discussions. We thank 
the anonymous referee for his/her very helpful suggestions which substantially improved the quality 
of this paper. We are supported in part by NSF AST-1109039, NNX13AE63G, NSFC 11173019, 
and the University of Missouri Research Board. 
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